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Abstract 

Metastability and long-range interactions of Mg, Ca, and Sr in the lowest-energy metastable 
3 P2 state are investigated. The calculated lifetimes are 38 minutes for Mg*, 118 minutes for 
Ca*, and 17 minutes for Sr*, supporting feasibility of cooling and trapping experiments. The 
quadrupole-quadrupole long-range interactions of two metastable atoms are evaluated for various 
molecular symmetries. Hund's case (c) A g potential possesses a large 100-1000 K potential barrier. 
Therefore magnetic trap losses can possibly be reduced using cold metastable atoms in a stretched 
M = 2 state. Calculations were performed in the framework of ab initio relativistic configuration 
interaction method coupled with the random-phase approximation. 



PACS numbers: 31.10,+z, 34.20. Cf, 32.10.Dk, 31.15.Ar 



I. INTRODUCTION 



This work is motivated by emerging experiments on cooling and trapping of alkaline-earth 
atoms ( see, e.g., Ref. [0], 0, |[]). In particular, the long-lived lowest-energy 3 P 2 state can 
serve as an effective ground state in such experiments. Recently Katori et al. J3J cooled and 
trapped metastable Sr and proposed a scheme for a continuous atom laser with a possible 
application in atom lithography. The purpose of this work is to evaluate properties of 
metastable nsnp 3 P 2 states for Mg (n = 3), Ca (n = 4), and Sr (n = 5). In particular, we 
calculated decay rates of the nsnp 3 P2 states. The resulting lifetimes are on the order of 
10 3 — 10 4 s supporting the feasibility of the experiments. 

Ultracold collision properties, including scattering lengths, are sensitive to long-range 
atomic interactions. The dominant van der Waals interaction of two atoms in their respective 
3 P 2 states is described in terms of the quadrupole moment of the atomic state. To assist in 
determining molecular potentials, the atomic quadrupole moments of Mg*, Ca*, and Sr* are 
also calculated here, and the relevant C5 coefficients for various molecular symmetries are 
tabulated. These coefficients are substantially larger compared to those for metastable noble- 
gases f|. In particular, Hund's case (c) 4 5 potential possesses a large 100-1000 K potential 
barrier. Therefore magnetic-trap losses can possibly be reduced using cold metastable atoms 
in a stretched M = 2 state || . 

II. METHOD 

Ab initio relativistic valence configuration-interaction (CI) method coupled with random- 
phase approximation (RPA) was employed here. A detailed description of this method |7j 
will be published elsewhere; only a brief discussion is presented here. In this method the 
wave-function is expanded in terms of two-particle basis functions as 

*(7rJM)=5>fci$«(7rJM) . (1) 
k>i 

Here J is the total angular momentum with projection M, and ir is the parity of the state 
The weights Cki and energies are found by solving the eigen- value problem based on the 
no-pair Hamiltonian ||. The basis functions are defined in the subspace of virtual orbitals 

( vrJM) = r, kl £ (fl,.,,,,, a\ nkKkmk} a\ niKimi} \0 core ) , (2) 

m k mi 

where sets {nnm} enumerate quantum numbers, 77^ = 1 — \5 nkni 5 KkKl is a normalization 
factor, at are creation operators, and the quasi- vacuum state |0 core ) corresponds to a closed- 
shell core. The one-particle orbitals are found in the-frozen core (V N ~ 2 ) approximation, 
i.e. the Dirac-Hartree-Fock (DHF) equations are solved self-consistently for core orbitals, 
and the virtual orbitals are determined using the resulting DHF potential of the core. The 
employed set of basis functions in Eq. @ is essentially complete, and one could interpret 
the solution of the eigen-value problem as treating the strong Coulomb repulsion of the two 
valence electrons to all orders of perturbation theory. This valence-CI method being exact 
for He-like systems 0, represents an approximation for the alkaline-earth systems. The 
core-polarization effects and the Breit interaction are neglected here. 

Once the wave-functions are determined, the matrix elements of a one-particle opera- 
tor Z = Y.v 

computed by forming products (^f\Z\^/j). The additional RPA 
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FIG. 1: Single-photon decay channels from the lowest-energy 3 P 2 state of alkaline-earth atoms. 



approximation [10|, describing shielding of an externally applied field by core electrons, con- 
stitutes a substitution of the "bare" one-particle matrix elements 2^ by "dressed" matrix 
elements zf^ A . Such an approach sums a certain class of many-body diagrams to all orders 
of perturbation theory. 

The configuration-interaction eigenvalue problem was solved numerically using the B- 

nn 



spline basis set technique 
h 



The employed basis set included the partial waves Si/ 2 - 
u/2- For each partial wave, 40 positive-energy basis functions approximated by B-splines 
represented a complete set. An inclusion of the 30 lowest-energy basis functions for each 
partial wave was found to be sufficient for the goals of this work. 



III. 



LIFETIME OF 3 P 2 STATES 



Single-photon transitions from the metastable 3 P 2 state are presented in Fig. |l|. It will 
be demonstrated that for Mg and Ca this level predominantly decays through magnetic- 
quadrupole (M2) transition to the ground state. For Sr the main decay channel is a 
magnetic-dipole (Ml) transition to the J = 1 level of the same 3 P fine-structure multiplet. 
Estimates show that the magnetic-octupole (M3) 3 P 2 — > 3 P± transition and all second-order 
two-photon decays contribute a negligibly small fraction to the total decay rate. Further, 
cold-atom experiments are focused on isotopes without nuclear spin to avoid negotiating 
molecular potentials complicated by hyperfme structure; we will not consider otherwise 
important electric-dipole hyperfine-structure induced decays [1^ . 

The data on magnetic-quadrupole transitions (M2) are scarce. For Mg the 3 P 2 — 1 So M2 
rates were estimated by Mizushima [13| and Garstang [l2| more than three decades ago and 



recently by Jonsson and Fischer [TJj]. There are no published data on M2 transition rates 
for Ca and Sr. 

A general treatment of multipole electromagnetic transitions in the relativistic frame- 
work can be found, for example, in Ref. ||. The Einstein coefficient Am 2 for a magnetic- 
quadrupole transition / — > F is given byfT9fl 
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A-M2 = ^-rCt 7 00 5 



|QMA4 (2) ||ttj) 
2J/ + 1 



(3) 



where u is the photon frequency and a is the fine- structure constant. The relevant one- 
particle reduced matrix element can be expressed as 
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r { 15 Wf) r ^ G ^ F ^ + F ^(r)]dr. (4) 

Here is the normalized spherical harmonic [15|], G(F) are large (small) radial components 
of a wavef unction, k = (I— j)(2j + l), and photon wavevector is defined as k = ui/c. The term 
in curly brackets approaches unity in the long-wavelength approximation. Although the full 
retarded form of the matrix element was used in the calculations, the long-wavelength limit 
is well satisfied for transitions between the lowest-energy valence states of alkaline-earth 
atoms. While the matrix element is largely independent of the transition frequency u, the 
rate A M2 depends on it strongly ( ~ u> 5 ) and the experimental energy intervals were used 
in the final tabulation. 

The calculated CI+RPA magnetic-quadrupole transition rates are presented in Table |. 
It was found that the RPA-dressing of matrix elements affects these rates only on the order 
of 1%. The calculated M2 transition rate for Mg*, 4.41 x 10~ 4 1/s, is in a fair agreement 



with the multiconfiguration Dirac-Fock result, 3.983 x 10 4 1/s, by Jonsson and Fischer [14 
Both rates are as twice as large compared to the earlier estimate [J2], 1.6 x 10~ 4 1/s. 



Magnetic-dipole 3 P 2 -+ 3 Pi transitions occur between the levels of the same fine-structure 
multiplet. The associated decay rate can be reliably estimated in the nonrelativistic approx- 
imation as 

A M1 = -a 5 uj 3 . (5) 
o 

The Ml rates, presented in Table |, were computed using experimental energy intervals. 
The determined rate for Mg is 3% larger than in calculations [Oj] where theoretical fine- 
structure splittings were employed. The discrepancy is accounted for by a difference between 
theoretical and experimental energies. 

We also calculated electric-quadrupole (E2) decay rates to the 3 Pi i0 fine-structure levels. 
The results listed in Table | show that the contribution of these decay channels is small 
compared to the M2 and Ml rates. 

From the analysis of Table [I] one finds that for Mg and Ca the 3 P 2 state predominantly 
decays through magnetic-quadrupole (M2) transition to the ground state. For Sr the 
main decay channel is a magnetic-dipole (Ml) transition to the J = 1 level of the same 3 P 
fine-structure multiplet. Overall, the calculated lifetimes for Mg*,Ca*, and Sr* are on the 
order of 10 3 — 10 4 seconds, thus favoring the usage of these metastable alkaline-earth atoms 
in cooling and trapping experiments. 



IV. QUADRUPOLE-QUADRUPOLE LONG-RANGE INTERACTIONS 

Ultracold collision properties, including scattering lengths, are sensitive to long-range 
atomic interactions. Here we focus on the long-range interactions of two alkaline-earths 
atoms in their respective 3 P 2 states. At large internuclear separations, R, atomic wavefunc- 
tions are perturbed by the axially-symmetric molecular potential, which itself depends on 
the electronic coordinates and the internuclear distance. This potential can be expanded in a 
sum of interactions of various atomic multipoles [IB|. The lowest-order correction to molec- 



ular term energies correlating to two 3 P 2 atomic states arises from a quadrupole-quadrupole 
contribution jT7| 

1 2 4' 

VQQ = # E a ( 2 - m )!(2 + a*)I {Q " )l iQ -" )n ' (6) 
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Here subscripts I and II denote the subspaces of the electronic coordinates of the two atoms, 
and the quadrupole spherical tensor is defined as 



Q^HelErfCfft), (7) 

i 

where the summation is over atomic electrons, is the position vector of electron i, and 
Cjpfii) are normalized spherical harmonics fll~5| . 

The quadrupole-quadrupole interactions are parameterized by the quadrupole moment 
Q of the 3 P 2 atomic state, defined conventionally as 

Q = 2 ( 3 P 2 ,Mj = 2\Q \ 3 P 2 ,Mj = 2). (8) 

This quadrupole moment is related to the reduced matrix element of the tensor, Eq. (|7|), as 
Q( 3 f > 2) = (^HQH 3 ^)- The associated one-particle reduced matrix element is given by 

(MQWh) = («il|C (2) lh-> / AG^G^r) + F t (r)F J (r)]dr. (9) 

<J 

Using the valence Cl-method coupled with the RPA dressing of matrix elements we cal- 
culated the quadrupole moments (see Table 0). It was found that the RPA sequence of 
diagrams reduces the final result only by 0.1% for Mg, 0.3% for Ca, and 0.5% for Sr. Due to 
the diffuse nature of valence states, these quadrupole moments are significantly larger than 
those found for metastable noble-gas atoms [El], where the hole in the outer p 3 / 2 subshell 
determines the quadrupole moment. 

As in the case of metastable 3 P 2 noble-gas atoms |5|, I8fl , the long-range molecular po- 
tentials of metastable alkaline-earth atoms depend on their spatial orientation, i.e. the 
interactions are anisotropic. Altogether, 15 distinct molecular states correlate to the two 
atomic 3 P 2 states at large separations. The first order correction to a molecular term may 
be represented as 

17«(J2) = §. (10) 



The relevant constants C5 for various molecular symmetries are given in Table |TJ, where 
the states are characterized using Hund's case (c). These coefficients were obtained by 
diagonalizing the quadrupole-quadrupole interaction, Eq. in the basis of products of 
atomic states \ 3 P 2 , M)i Cg> \ 3 P 2l M')n for a given symmetry Q satisfying M + M' = Q. 

Seven of the resulting quadrupole-quadrupole long-range potentials are attractive at large 
distances and eight are repulsive. We present the most repulsive Ojj" and the most attractive 
3 g potentials in Fig. [2[ It is worth discussing the repulsive 4 9 potential also shown in Fig. |2|. 
Provided the metastable atoms in a magnetic trap are prepared in a stretched M = 2 state, 
the collisions would occur along this repulsive VI — 4 potential. Although the medium-range 
part of the potential is attractive, the resulting barrier will be on the order of 100-1000 K 
high. Therefore ultracold collisions of the 3 P 2 alkaline-earth atoms in a stretched M = 2 
state can be effectively shielded from the losses occurring at small internuclear distances ||. 
This barrier is more pronounced compared to metastable noble-gas atoms ||, where the 
potential barriers were a few nK high. Clearly, calculations of the second-order C§ dispersion 
coefficients and intermediate- and short-range parts of the potentials are needed for a more 
quantitative description. 
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TABLE I: Einstein coefficients A for decays from the lowest-energy 3 i"2 states in s 1 . Notation x[y] 
stands for x x 10 y . 



Type, final state Mg Ca Sr 

M2,ns 2 % 4.41[-4] 1.25 [-4] 1.27[-4] 

Ml, nsnp^Pi 9.12[-7] 1.60[-5] 8.26[-4] 
E2, nsnp 3 Pt 1[-12] 3[-10] 3[-7] 
E2, nsnp 3 P 3 [-12] l[-9] l[-6] 

A, total 4.42[-4] 1.41 [-4] 9.55[-4] 



TABLE II: Quadrupole moments Q of the lowest-energy metastable 3 P2 states in a.u. 

"Mg Ca Sr 
8.59 13.6 16.4 



V. CONCLUSION 

To address the needs of emerging experiments on cooling and trapping of alkaline-earth 
atoms, we performed ah initio relativistic calculations of lifetimes and quadrupole moments 
of metastable Mg, Ca, and Sr. The determined lifetimes are 38 minutes for Mg*, 118 minutes 
for Ca*, and 17 minutes for Sr*, supporting the feasibility of experiments. In addition, we 
investigated long-range quadrupole-quadrupole interactions for molecular potentials corre- 
lating to two metastable atoms. Several resulting potentials possess pronounced barriers, 
which could be exploited to minimize trap losses. 

Thanks are due to H. Katori, C. Oates, and F. Riehle for stimulating discussions and V. 
Davis for careful reading of the manuscript. The developed numerical code was partially 
based on programs by Notre Dame group led by W.R. Johnson. This work was supported in 
part by the Chemical Sciences, Geosciences and Biosciences Division of the Office of Basic 
Energy Sciences, Office of Science, U.S. Department of Energy. 



TABLE III: C5 coefficients in units of Q 2 for molecular states ( Hund's case (c)) asymptotically 
correlating to two 3 i-2 state atoms. The quadrupole moments Q are listed in Table |l], and long- 
range interaction potentials are parameterized as U(R) = C5/-R 5 . 
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FIG. 2: The most repulsive 0+ and the most attractive 3 g long-range quadrupole-quadrupole 
interaction potentials correlating to two metastable 5s5p 3 P2 Sr atoms. The 4 g potential is also 
shown. 1 K = 0.69503 cm -1 . 
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